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The proton spin dependent structure function g\ is analysed using the unified 
. scheme incorporating both Altarelli-Parisi evolution and the double ln 2 (l/x) 

effects at low x. The latter are found to be very important in the region of 
x which can possibly be probed at HERA. Predictions for the polarized gluon 
distribution AG(x) at low x are also given. 
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It has recently been pointed out that the spin dependent structure function g\ at low x 
is dominated by the double logarithmic ln 2 (l/x) contributions i.e. by those terms of the 
perturbative expansion which correspond to the powers of ln 2 {l/x) at each order of the 
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expansion [|. Those contributions go beyond the LO or NLO order QCD evolution of 
polarised parton densities and in order to take them into account one has to include 
the resummed double ln 2 (l/x) terms in the coefficient and splitting functions [|J. In 
ref. U we have developed an alternative approach based on unitegrated distributions 
which for simplicity was formulated for the non-singlet distributions. In this note we 
wish to present the preliminary results for the singlet structure functions concentrat- 
ing on the proton polarized structure function gi(x,Q 2 ). Our main aim is to study 
gi(x, Q 2 ) in the region of the small values of x which can possibly be probed at HERA 
H and confront obtained results with those obtained from the LO Altarelli-Parisi evo- 
lution. We will also show similar results for the polarised gluon distribution AG(x, Q 2 ). 



The dominant contribution to the double ln 2 (l/x) resummation comes from the 

ladder diagrams with quark and gluon exchanges along the ladder. In what follows we 
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shall neglect for simplicity possible non-ladder bremsstrahlung terms which are rela- 
tively unimportant [fll 0. 



It is convenient to introduce the unintegrated (spin dependent) parton distributions 
fi(x',k 2 ) (i = u v ,d v ,u,d,s,g) where k 2 is the transverse momentum squared of the 
parton i and x' the longitudinal momentum fraction of the parent proton carried by 
a parton. The conventional (integrated) distributions Api(x,Q 2 ) are related in the 
following way to the unintegrated distributions fi(x', k 2 ): 

,„, rW 2 Ah? k 2 

A Pl (x,Q 2 ) = Ap { ?\x)+j k2 ^fi(x' = x(l + ^- 2 ),k 2 ) (1) 

where Ap| (x) is the nonperturbative part of the of the distributions and 

W 2 = Q 2 (- - 1) (2) 

x 

where as usual x is the Bjorken parameter and Q 2 = —q 2 with q denoting the four 
momentum transfer between leptons. The parameter k% is the infrared cut-off which 
will be set equal to 1 GeV 2 . The origin of the nonperturbative part Ap\ (x) can be 
viewed upon as originating from the non-perturbative region k 2 < k 2 ,, i.e. 

Apf\x)= J^f^k 2 ) (3) 

The spin dependent structure function gi(x, Q 2 ) of the proton is related in a standard 
way to the (integrated) parton distributions: 



\{Au v {x : Q 2 ) + 2Au{x, Q 2 )) + \{Ad v (x, Q 2 ) + 2Au{x, Q 2 ) + 2As(x, Q 2 )) 
9 9 



(4) 

where Au v (x, Q 2 ) = Ap Uv (x,Q 2 ) etc. We assume Am = Ad and confine ourselves to 
the number of flavours Np equal to three. 



The valence quarks distributions and asymmetric part of the sea 

fus{x\k 2 ) = U{x\k 2 )-f- s {x',k 2 ) (5) 

will correspond to ladder diagrams with quark exchange along the ladder. The singlet 
distributions 

f s (x', k 2 ) = f Uv (x', k 2 ) + f dv (x', k 2 ) + 4f a (x', k 2 ) + 2f- s {x' , k 2 ) (6) 



and the gluon distributions f g (x',k 2 ) will correspond to ladder diagrams with both 
quark (antiquark) and gluon exchanges along the ladder. 



The sum of double logarithmic ln 2 (l/x) terms corresponding to ladder diagrams is 
generated by the following integral equations: 



f k (x\k 2 ) = tf\x',k 2 ) + ^AP qq (0) f 

Z7T Jx> Z Jk 

(k = u v , d v , US) 



^dz r^i-dk' 2 fx' , 



(7) 



fs(x',k 2 )=jf(x',k 2 )+ 



a q r^dz r k2 / z dk' 2 



2n Jx> z Jkl 



k' 2 



X ' "2 \ , ad fn\ t ( X 7J2 



&P q Mfs ~, k' 2 + AP gg (0)f g -, k 



f g (x',k 2 ) = ff\x',k 2 ) + 



a. 



^dz r k2 / z dk' 2 



k' 2 



2n Jx' z Jkl 



x' 



^P gq (0)fs{j,k' 2 j+AP gg (0)f g ^-,k' 2 ^ 

with APjj(O) = APij(z = 0) where AP^z) denote the LO splitting functions describing 
evolution of spin dependent parton densities and the inhomogeneous terms f^(x,' , k 2 ) 
will be specified later. To be precise we have: 

AP w (0) = \ 
AP qg (0) = -N F 

AP gg (0) = 12 (9) 

Equation (|7D is similar to the corresponding equation in QED describing the dou- 
ble logarithmic resummation generated by ladder diagrams with fermion exchange |7| . 
The problem of double logarithmic asymptotics in QCD in the non-singlet channels 
was also discussed in ref ||. 

Equations (|7|, generate singular power behaviour of the spin dependent parton 
distributions and of the spin dependent structure functions g\ at small x i.e. 



g? s (x,Q 2 
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g s Ax,Q 2 ) 

AG(x,Q 2 ) ~ x~ Xs (10) 

where g± s = g\—g\ and gf = gi+g™ respectively and AG is the spin dependent gluon 
distribution. This behaviour reflects similar small x' behaviour of the unintegrated 
distributions. Exponents Xjsrs S are given by the following formulas: 



A 



NS 



A, 



2tt 



AP w (0) 
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where 
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AP qq (0) + A? M (0) + J(AF„(0) - AP gg (0)) 2 + 4AP g9 (0)AP OT (0) 



(12) 



Both equations (|7|||) and the exponents \ns,s acquire additional contributions from the 
non-ladder diagrams which are however relatively small |], gj. The power-like behaviour 
(O) remains the leading small x behaviour of the structure functions provided that 
their non-perturbative parts are less singular. This takes place if the latter are assumed 
to have the Regge pole like behaviour with the corresponding intercept (s) being near 0. 



Following ref. || we extend equations ([j],[8]) and add to their right hand side the 
contributions coming from the remaining parts of the splitting functions APy(^). We 
also allow the coupling a s to run setting k 2 as the relevant scale. In this way we 
obtain unified system of equations which contain both the complete LO Altarelli-Parisi 
evolution and the double logarithmic ln 2 (l/x) effects at low x. The corresponding 
system of equations reads: 



f k (x',k 2 ) = fl°\x',k 2 ) + 



a s (k 



2 M r l dz 
2n 3 Jx' z 



k 2 /*dk' 2 fx' " 



a s (k 2 ) r^dk' 2 A r±dz(z + z 2 )f k (^,k' 2 )-2zf k (x',k' 2 ) 



2tt 



(yk — w^j, U S) , 



k' 2 3Jx> z 
a s (k 2 ) r* 2 dk' 2 



2tt 



k' 2 



2 + -ln(l - x' 



1-z 

f k (x',k' 2 ) 



(13) 



fs{x\k 2 ) = ff\x\k 2 ) + 



a s (k 2 ) r±dz r^hdk' 2 A 



2tt 



x 1 Z 



X 



l- 



+ 



(k 2 ) r k2 dk' 2 A /-i dz(z + z 2 )f s (^,k' 2 )-2zf s (x',k' 2 ) 



2tt Jul k' 2 3J X > z 
a s (k 2 ) r* 2 dk' 2 



2vr 



2tt Jul k' 2 
^dz r k2 /* dk' 2 



2 + -ln(l -x 



x > z h\ k' 2 



u 



x 



J2 



1 - Z 

fs(x',k' 2 )+ 
fc2 dk' 2 r 1 dz 



/,•;; k' 2 J,' : 



.r 



2zf g (-k r2 ) 



f g (x',k 2 ) = f^(x',k 2 ) + 



a. 



(A; 2 ) /-i^ r^/zdk' 2 



2tt 



* y fe 2 fc' 2 3 



2vr 



k2 dk' 2 r^dz A, 
kl "k 12 Jx> T^~3' 

(k 2 ) f dk' 2 fldz 



idz r k2 /zdk' 2 



z hi k' 2 



X 



2vr Jkl k' 2 Jx 
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f 9 (^k' 2 )-f g (x',k' 2 ) ^ 



a s {k 2 ) f k2 dk' 2 



2tt M k' 2 



1-z 



y-y + 6Ml-^') 



fg(X } k ) 



(14) 



The inhomogeneous terms f^°\x',k 2 ) are expressed in terms of the input (inte- 
grated) parton distributions and are the same as in the case of the LO Altarelli Parisi 
evolution [O: 



fl°\x',k 2 ) 



a s (k 2 ) 4 /•! dz (1 + z 2 )Api 0) (^) - 2zA P f{x' 



Mi 



2n 3 A' z 
a s (k 2 



2ty 



2 + -ln(l-x') 



1 - z 
Apf{x') 



(15) 



a s (k 2 



2ti 3 A' z 



a s {k 2 ) 4 /•! ^ (1 + z 2 )Apf\f) - 2^A4 0) (a; / 



,(°) ( 



2vr 



1 - ^ 

1 dz 



(2 + -Jn(l - x'))Apf(x') +N F f — (1 - 2z)Apf 
3 jx' z 



a s (A; 2 ) 



4 ^dz 
3 
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2tt 

a s (/c 2 ) r 1 (is 



/ -(2 - ^)A^ 0) (-) + - ^ + 6/n(l - x'))Apf 



2tt 



6 — 



1 — 2 y z 



(16) 



Equations ( |T3l , |T4D together with C[T5| , p~6|) and (0) reduce to the LO Altarelli-Parisi 



evolution equations with the starting (integrated) distributions Ap®(x) after we set the 
upper integration limit over dk' 2 equal to k 2 in all terms in equations (pT3|JT^ ) and if we 
set Q 2 in place of W 2 as the upper integration limit in the integral in eq. (fj). 

We solved equations ( |T3| , [T4|) assuming the following simple parametrisation of the 
input distributions: 

Apf (x) = Ni(l - x) v > (17) 

where we set r] Uv = 77^ = 3, r] a = 7] s = 7 and r) g = 5. The normalisation constants 
Ni were determined by imposing the Bjorken sum-rule for Am^ — Ad^ and requiring 
that the first moments of all other distributions are the same as those determined from 
the recent QCD analysis 0. All distributions Apf\x) behave as x° in the limit x — > 
that corresponds to the implicit assumption that the Regge poles which correspond to 
axial vector mesons, which should control the small x behaviour of g\ [ffO] have their 
intercept equal to 0. We checked that the parametrisation fllTD combined with equa- 
tions (ffl|^Jl3lJl4]) gives reasonable description of the recent SMC data on gi(x, Q 2 ) [l ll . 



In Fig.l we show gi(x, Q 2 ) for Q 2 = lOGeV 2 in the small x region which can pos- 



sibly be probed at HERA. We show predictions based on the equations (|l3l , |i4 |l|) and 
confront them with the expectations which follow from solving the LO Altarelli-Parisi 
evolution equations with the input distributions at Ql = lGeV 2 given by equation(s) 
(|17|). We also show in this Figure the "experimental" points which were obtained from 
the extrapolations based on the NLO QCD analysis together with estimated statisti- 
cal errors 0. We see that the structure function gi(x,Q 2 ) which contains effects of 
the double ln 2 (l/x) resummation begins to differ from that calculated within the LO 
Altarelli Parisi equations already for x ~ 10~ 3 . It is however comparable to the struc- 
ture function obtained from the NLO analysis for x > which is indicated by the 
"experimental" points. This is presumably partially an artifact of the difference in the 
input distributions but it also reflects the fact that the NLO approximation contains 
the first two terms of the double ln 2 (l/x) resummation in the corresponding splitting 
and coefficient functions. It can also be seen from Fig.l that the (complete) double 
ln 2 (l/x) resummation generates the structure function which is significantly steeper 
than that obtained from the NLO QCD analysis and the difference between those two 
extrapolations becomes significant for x < 10~ 4 . 



In Fig. 2 we show our predictions for the polarized gluon distribution AG(x, Q 2 ) = 
Ap g (x,Q 2 ) for Q 2 = lOGeV 2 and confront it with the polarized gluon distribution 
obtained from the LO Altarelli-Parisi equations. The gluon distributions exhibits char- 
acteristic x~ Xs behaviour with \$ ~ 1. Similar behaviour is exhibited by the structure 
function gf(x, Q 2 ) itself. 

To sum up we have presented results of the analysis of the "unified" equations 
which contain the LO Altarelli Parisi evolution and the double ln 2 (l/x) effects at low 
x. As the first approximation we considered those double ln 2 (l/x) effects which are 
generated by ladder diagrams. The double logarithmic effects were found to be very 
important and they should in principle be visible in possible HERA measurements (cf. 
Fig. 1). 
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Figure Captions 

Fig.l The structure function cfi(x, Q 2 ) for Q 2 = lOGeV 2 plotted as the function of x. 
Solid line represents this structure function with the double ln 2 (l/x) terms included 
and the dashed line corresponds to g\ obtained from the LO Altarelli-Parisi equations 
The "experimental" points are based on the NLO QCD predictions with the statistical 
errors expected at HERA 

Fig. 2 The spin dependent gluon distribution AG(x, Q 2 ) for Q 2 = lOGeV 2 plotted 
as the function of x. Solid line represents AG(x, Q 2 ) with the double ln 2 (l/x) terms 
included and the dashed line corresponds to the AG(x, Q 2 ) obtained from the LO 
Altarelli-Parisi equations. 
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